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Abstract 

In this paper we study three different functional approaches to classical thermal field theory, 
which turn out to be the classical counterparts of three well-known different formulations of quan- 
tum thermal field theory: the Closed-Time Path (CTP) formalism, the Thermofield Dynamics 
(TFD) and the Matsubara approach. 
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I. INTRODUCTION 



In recent years, there has been a growing interest in the physics of collisions between 
heavy ions [1, 2]. In particular, one of the most interesting topics in the field has become the 
study of the quark-gluon plasma (QGP), which is a state of baryonic matter characterized 
by high values of temperature and density and by the fact that quarks are not confined [3]. 
A better understanding of the physics of this state would be an important step forward in 
high energy physics and it could have important implications for early universe cosmology 
as well [4]. 

The theoretical framework which is required to study the QGP is the finite temperature 
quantum field theory. In fact, this formalism can take into account the quantum, the 
relativistic and the statistical features of the system at the same time. There are evidences 
that in the QGP gluons are characterized by very large occupation numbers. In this state, 
one would expect the system to be described by a classical field theory [5]. In fact, it is 
known that QED is well-approximated by the classical Maxwell theory when occupation 
numbers of photons become large. Analogously, when occupation numbers of gluons are 
large, it should be possible to approximate QCD with a classical gauge field theory [6]. 

These considerations have led various authors to study classical field theories at high 
temperatures and, to our knowledge, three main approaches have been pursued in the lit- 
erature so far. The first approach is due to Aarts and collaborators [7-11], who computed 
classical thermal correlation functions {(f>(x) ■ ■ ■ <f>(x'))p by solving the Hamiltonian equations 
of motion and averaging the solutions over a canonical ensemble of initial conditions. An 
alternative approach has been proposed by F. Cooper et al. [12] and it is based on the op- 
eratorial formalism developed by P.C. Martin, E.D. Siggia and H.A. Rose [13] in the 1970s 
and known as MSR formalism. Finally, S.Jeon [14] has recently studied a classical scalar 
field theory at finite temperature by resorting to the path integral approach to classical 
mechanics proposed in ref. [15] at the end of the 80s. 

Despite the formal differences, there is a common feature shared by the three classical 
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approaches mentioned above: they are all related to the quantum Closed Time Path (CTP) 
formalism [16, 17] in the high temperature limit. However, quantum field theory at finite 
temperature admits several other equivalent formulations, the most popular ones being the 
Matsubara formalism [18] and the Thermofield Dynamics (TFD) [19]. The goal of this paper 
is to study the "classical analogs" of the TFD and the Matsubara formalism as well. To 
this end, we will present a general framework for classical thermal field theory based on 
the operatorial approach to classical mechanics proposed in the 30 's by Koopman and von 
Neumann (KvN) [20-22] and on the associated classical path-integral (CPI) formulation 
[15]. We believe that these tools (KvN and CPI) are the most suitable ones for studying the 
interplay between classical and quantum mechanics [23] either with or without temperature. 
Within this theoretical framework, we will show how the classical finite temperature field 
theory can be formulated in three different ways, which are in a certain sense the " classical 
counterparts" of the three quantum formalisms mentioned above: the CTP, the TFD and the 
Matsubara one. All these formalisms will be illustrated by considering the simple example 
of a scalar field with quartic self-interaction. 

The paper is organized as follows: after reviewing briefly the various equivalent ap- 
proaches to quantum thermal field theory (Section II) and the KvN and CPI formulations 
of classical mechanics (Section III), in Section IV we will present the approach to classical 
thermal field theory developed in [14]. From a physical point of view, this is the most in- 
tuitive approach to classical thermal field theory, as it is based on the same principles used 
by Aarts and Smit in [7, 8], i.e. on solving the equations of motion and averaging over the 
initial conditions using the canonical distribution. As we already mentioned, this approach 
turns out to be related to the CTP approach to quantum thermal field theory in the high 
temperature limit. 

Then, in the second part of the paper we shall develop two new approaches to classical 
thermal field theory which display both a remarkable formal analogy and a quantitative 
agreement with their quantum counterparts at high temperatures, i.e. the TFD approach 
and the Matsubara formalism. In particular, in Section V we will study the TFD approach 
to classical thermal field theory, which is based on the idea of implementing thermal aver- 
ages as expectation values on a particular state in the classical Hilbert space of ref. 
[24]. In Section VI we will study the classical counterpart of the quantum Matsubara for- 
malism. Both quantum and classical Matsubara approaches allow us to compute only static 



4 



properties of system, since the time variable is formally restricted to the purely imaginary 
axis. However, while in the quantum case time is an element of a complex plane, in the 
classical case time is an element of a complex superspace [25]. As we will see, the appearance 
of superspace is just a natural consequence of the fact that the path integral formulation 
of classical mechanics features a universal N = 2 supersymmetry and, as such, admits a 
very compact representation in terms of some suitably defined superfields [23]. Finally, 
some further background material as well as some technical details were included in two 
appendices. 



II. QUANTUM THERMAL FIELD THEORY 

In order to make the paper as self-contained as possible, in this section we will review 
some aspects of quantum field theory at finite temperature. For a more exhaustive treatment 
we refer the reader to the review articles [26, 27] or the textbooks [28, 29]. The main goal of 
this section is to show how perturbative calculations can be implemented in several different 
ways, and to derive the corresponding sets of Feynman rules. For concreteness, the main 
concepts will be illustrated by considering the simple case of a scalar field described by the 
Hamiltonian 

H = J* x $ + Zf + ^ + i£}, (1) 

where n is the canonical momentum associated with 0. 

In the case of a scalar field in thermal equilibrium at temperature T = 1/(3, the quantities 
of interest are the thermal correlation functions, which are defined as expectation values of 
time-ordered products of fields 0: 



Trje-^T^!)---^)]} 



G fi (x!, ...,x n ) = (T[(f)(xi) ■ ■ ■ 4>{x n )\)p = Tre-PH ' ^ 

where T[- ■ •] stands for the time-ordered product and p = e~^ H is the statistical operator 
describing a canonical ensemble. These quantities are clearly a generalization of the corre- 
lation functions one encounters in ordinary quantum field theory, since they reduce to the 
vacuum expectation value of a time-ordered product of fields in the zero-temperature limit 

/3 ->■ oo: 

G/ 3 {x 1 ,...,x n ) = ^ — — > {0\T[<f>{xi)---<f>{x n )]\0). (3) 



Before deriving a set of Feynman rules which will allow us to calculate thermal correlation 
functions in a systematic way, it is useful to derive an important result known as Kubo- 
Martin-Schwinger (KMS) relation [30, 31]. This result is based on the simple idea that, if one 
allows time to be a complex variable, then the statistical operator e~^ H can be interpreted 
as a time evolution operator with a purely imaginary argument —i/3. As a consequence, the 
thermal average of the product of any two observables A and B in the Heisenberg picture 
satisfies the following identity: 

(AmV))B = jiiL^mm = ^mMt+m = W)Mt + (4) 

\_re~ ptl rce~ ptl 
This identity is the celebrated KMS relation. We should point out that this relation was 
obtained by using the cyclicity of the trace, and this property may fail if the trace is not 
finite [28]. We will encounter an example of such a behavior in Section VI. 

A. Path Integral Formulation 

The idea of interpreting the statistical operator e~ pH as a time evolution operator not 
only allow us to derive an important identity such as the KMS relation, but it also leads 
naturally to a path integral expression for the canonical partition function Zp = Tre^ 1311 
[32, 33]. In fact, by expanding the trace on the basis of the eigenstates of the field operator 
0(t, x) at a give time t , we get 

Z p = y^](<M o |e^|0,t o > = y^](0,* o -^|<Mo>, (5) 

where [d(f>] denotes a path integral over all the field configurations 0(to,x) with fixed time 
to- By resorting to the standard slicing procedure, the transition amplitude on the RHS of 
equation (5) can be easily expressed in terms of a path integral over all the field configurations 
that satisfy the periodic boundary condition 0(t o ,x) = 0(t o — i/?,x). In the particular case 
of a scalar field described by the Hamiltonian (1), we get 

Zp = Jv^e 13 , S = JdtJ^x^d^-^-^Y (6) 

where C is a path in the complex plane of time which starts at t and ends at t — ij3. 
Notice that the time to is arbitrary, and all that matters is that the initial and final points 
of the time contour C are separated by an amount —ij3. The path integral expression (6) 
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immediately suggests that thermal correlation functions evaluated at times along the path 
C could be systematically calculated starting from the following generating functional: 

/ V c <f) exp I iS + i / dt / d 3 x j<f) I 
Z[j] = J - ^ J -^-L i (7) 

V c (f) exp {iS} 

However, consistency requires that the time contour C lies within the domain of analyticity 
of the thermal correlation functions Gp(x\,--- ,x n ). This constraint is satisfied provided 
the time path C is nowhere directed upward in the complex plane of time [26]. Notice also 
that the time ordered product that appears in equation (2) refers now to the order along 
the time path C. 

Starting from equation (7), a set of Feynman rules can be derived as usual by treating the 
quartic coupling in (6) as a perturbation. In particular, the propagator can be determined by 
considering just the quadratic terms in equation (6). In this case, the generating functional 
can be calculated exactly and takes the form 



Zo\j] = exp j-^ J d 4 xd A x' j(x)Gp(x - x')j(a/)| , 



(8) 



where Gp(x — x') is the propagator, i.e. the thermal correlation function of two fields 0, 
which satisfies the following differential equation: 

(n c + m 2 )Gp(x-x') = -i5 c (t-t')5(x-x'). (9) 

Notice that, in this equation, both the differential operator D c and the Dirac delta 5 c (t — t') 
are defined along the time contour C. According to the definition (2) of thermal correlation 
functions, the propagator G^x — x') can be written as 

Gp(x - x) =G>(x- x') 9 c (t - t') + G < (x- x') e c (t' - t), (10) 

where G > (x — x') = (4>(x)(j)(x')) p, G^x — x') = (cf)(x')(f)(x))p and the Heaviside step function 
along the contour C is related to the Dirac delta in equation (9) by -^.0 c (t — t') — S c (t — t'). 
The two functions G > and G < are not independent of each other, since the KMS relation 
(4) implies that 

G > (t-t / ,x-x') = G<(t-t' + i/3,x-x'). (11) 
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The propagator Gp(x — x') can then be determined univocally by solving the differential 
equation (9) together with the boundary condition (11). After performing a Fourier transfor- 
mation with respect to the spatial variables only, the propagator in the (t, p) space reads [26] 

. . = G p (t - t', p) = \(e iE ^-^ + e f>E Ve -iE p (t-t>)\ Qc{t _ t>) + 

2Ep L V / 

+ (e-iE p (t-t>) + e PE Pe iE p (t-t>)^ _ t) j (12) 

where we defined E 2 = p 2 + m 2 , while n,B(E p ) = (e^p — l)" 1 is the Bose-Einstein density 
of states. The effect of the quartic coupling can then be taken into account as usual by 
introducing a vertex with the following rule in the (t, p) space: 



X ■< 



dt. (13) 



Although the Feynman rules (12) and (13) apply to any path C which is nowhere directed 
upward in the complex plane of time, concrete applications usually become simpler once a 
specific contour is chosen. In the remaining part of this section we will briefly review some 
of the most popular choices. 



B. Matsubara Formalism 



The first approach to quantum field theory at finite temperature was originally developed 
by Matsubara [18] in the 50's and it is based on a path the runs along the imaginary axis 
from t — to t — —i/3, as shown in Figure 1. Since we are ultimately interested in calculating 
thermal correlation functions at real times and the only intersection between the time contour 
in Figure 1 and the real axis is the starting point t — 0, this formalism will only allow us to 
calculate time-independent correlation functions of fields 0(0, x). As such, the Matsubara 
formalism is a good tool for studying static properties of a system at equilibrium, but it is 
not suited to the study of time-dependent phenomena. 

The path in Figure 1 can be conveniently parametrized by introducing the Euclidean 
time such that t = —He- The main advantage of working with the Euclidean time as 
opposed to the time variable t is that now belongs to the finite interval [0,13], and since 
a trivial consequence of the KMS relation is that Gp(— i/3, p) = Gp(0,p), the propagator 
Gp(—itE, p) is a periodic function of on such interval. For this reason, it can be expressed 
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as a Fourier series, and the Fourier coefficients can be easily calculated starting from the 
general expression (12) for the propagator. The result one gets is 

dt E G p {-it E ,v) e WntE = — 

'o 



G fi (u n ,p)= / dt E Gp(-it E ,p)e^ tE 
Jo 



(14) 



where n = 0, ±1, ±2, ... and the discrete frequencies u n = nir/f3 are known as Matsubara 
frequencies. Furthermore, in the (u n , p) space, the rule (13) for the vertex reduces to 




-9- 



(15) 



Clearly, one advantage of the Matsubara formalism is that its Feynman rules are fairly 
simple. In fact, for the simple model we are considering, they consist of only one propagator 
(see eq. (14)) and one vertex (see eq. (15)). Notice also that these Feynman rules bear a 
close similarity to the ones encountered in ordinary quantum field theory, the only difference 
being that the discrete Matsubara frequencies have now replaced the continuous component 
p° of the four-momentum. However, as we already mentioned, a severe limitation of this 
formalism is that it is not suitable for calculating dynamical properties of a system. 



Im(t) 



Re(t) 



FIG. 1: Purely imaginary time path associated with the Matsubara formalism. 
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FIG. 2: Time path associated with real time formalisms. 
C. Real Time Formalisms: CTP and Thermofield Dynamics 

The main reason why the Matsubara formalism does not lend itself to the calculation of 
time-dependent properties of a system is that the time contour in Figure 1 does not lie along 
the real axis. In fact, dynamical properties can be calculated much more easily by working 
with the time path C shown in Figure 2, which defines the so called real time formalisms. 
Such path starts at t — — T, ends at t — — T — i(3, and consists of four segments C{. The 
distance between the two horizontal segments C\ and C 2 is controlled by the parameter a. 
In the limit where T — >■ oo, the two vertical segments C 3 and C 4 are pushed to infinity and 
therefore we expect their contribution to become negligible. This statement can be made 
more quantitative by breaking down the propagator Gp into several different components 
according to the following definition: 

Gij(t - 1', P ) = G p (t - 1', P ) e - tlRe{t - t,)l with t e d, t' e C r (16) 

So, for example, G 13 stands for the propagator that "connects" times on the segment C\ 
with times on the segment C3, and so on. Notice that, in equation (16), we have introduced 
an exponential term depending on the positive, infinitesimal parameter e. The introduction 
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of this term is necessary to render the limit T — > oo well-behaved, and it is analogous to 
the "e convention" usually adopted in ordinary quantum field theory By combining the 
general expression (12) for the propagator Gp together with equation (16), it is easy to show 
that the propagators connecting times along horizontal segments with times along vertical 
segments vanish in the limit T — > oo, i.e. 

G13, G23, Gu, G24 — > 0. (17) 

Thus, horizontal and vertical segments effectively decouple in this limit and, since we are 
ultimately interested in calculating thermal correlation functions at times along the real axis 
(i.e. along the segment Ci), the vertical segments C 3 and C 4 can be neglected. 

At the same time, in the limit T — >• 00 the horizontal segments C\ and C 2 essentially 
reduce to "two copies" of the real axis which can be respectively parametrized as t 1 (s) = s 
and £2(3) = — s — icr(3, with —00 < s < 00. It is therefore possible to consider the Fourier 
transform of the propagators Gu, Gi 2 , G 2 \ and G 22 with respect to the parameter s, i.e. 

/oo 
dse-^'GvMslp), i,j = 1,2. (18) 

-00 

Notice that these Fourier transforms are well defined because of the exponential term in 
equation (16). By combining equations (12), (16) and (18) it is easy to derive the following 
rules for the propagators in the space of momenta 



• = Gn(p) = ~ % -^— + 2nn B (\p°\W - m 2 ) (19a) 

p 2 — m 2 + it 

=Gia(p) = 27re^° [0( -p°) + n B (\p°\)} 5(p 2 - m 2 ) (19b) 
. = G 21 (p) = 27re-^° [6(p°) + n B (\p°\)] 5(p 2 - m 2 ) (19c) 



^ . =G 22 (p) = — - + 2Tm B (\p°\)5(p 2 -m 2 ) (19d) 

After "splitting" the propagator Gp into the various components defined in equation (16), 
we can also split the integral along C appearing in the rule (13) for the vertex into a sum of 
integrals along the four segments C«. However, in the limit T — > 00 only the contributions 
coming from the segments C\ and C 2 are relevant, which leads to the following rules for the 
vertices: 

-ig V =ig. (20) 
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The sign difference between the two vertices is due to the fact that segments C\ and C 2 have 
opposite directions, as clearly shown in Figure 2. We see therefore that the Feynman rules 
for real time formalisms are more complicated than the ones for the Matsubara formalism 
since, for example, four propagators and two vertices are now needed to describe a scalar 
field with a quartic interaction. This effectively amounts to a doubling of the degrees of 
freedom necessary to describe the system, and it is a price we need to pay in order to study 
dynamical properties. 

In conclusion, we should mention that two popular choices for the parameter a appearing 
in the Feynman rules (19a) - (19d) for the propagators are a = and a — 1/2. In the 
first case, the segments C\ and C 2 form essentially a closed contour, which explains why 
the associated formalism is known as Closed Time Path (CTP) formalism [28]. The choice 
a = 1/2 reproduces instead the Feynman rules of an operatorial approach to quantum 
thermal field theory known as Thermofield Dynamics (TFD) [19]. Notice that, in this case, 
the Feynman rules for the propagators become slightly simpler, since it can be easily shown 
[28] that G 12 (p)=G 21 (p). 

III. HILBERT SPACE AND PATH INTEGRAL APPROACHES TO CLASSICAL 
MECHANICS 

In this section we will briefly review the Hilbert space approach to classical statistical 
mechanics, which was originally developed by Koopman [20] and von Neumann [21, 22] 
in the 1930s. Then, we shall also present its path integral counterpart [15]. In order to 
review the main results needed in this paper, we shall restrict ourselves to the simple case 
of a one-dimensional system, since the extension to systems with more degrees of freedom 
is straightforward. We refer the reader to the original literature for further details, and in 
particular to [34] for an application to classical Yang-Mills theories. 

A. Hilbert space approach of Koopman and von Neumann (KvN) 

The formalism developed by Koopman and von Neumann (KvN) is based on a very 
simple idea: instead of describing a classical statistical ensemble by means of a probability 
distribution p(<p,t) in phase space, where ip = (q,p) denotes the phase space variables, one 
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can use a complex function ip(ip,t) such that the probability distribution is given by the 
square of its modulus, i.e. 

P (<p,t) = m^t)\ 2 . (21) 

Since the probability distribution is normalized, ip(tp, t) must be a square integrable function 
and therefore it belongs to a Hilbert space, which will be called KvN space. The functions 
i/)((p,t) will be called instead KvN waves in order to distinguish them from the quantum 
wave functions. Switching to the Dirac notation, tp(ip,t) can be rewritten as (ip\ip,t). Note 
in fact that the eigenstates of the phase-space operators \ip) = \q,p) form a basis of the KvN 
space. This is due to the fact that, in classical mechanics, the operators q and p commute, 
i.e. 

[<p a ,<p b \ = 0, (22) 

and therefore they can be simultaneously diagonalized. 

There is however a price to pay in passing from a description of classical ensembles in 
terms of probability distributions to a description involving states in a Hilbert space. In 
fact, according to equation (21), every distribution p(<p,t) can be described by a KvN wave 
of the form: 

M<P,t) = y/faj)e*to*\ (23) 

where ee((p, t) is an arbitrary local phase. By varying the local phase we get an infinite set of 
distinct states in the KvN space which are all associated with the same distribution in phase 
space, and therefore have the same physical content. Thus, unlike in quantum mechanics, 
the local phase of a KvN state does not carry any physical information. In fact, it was 
shown in [35] that the phase and the modulus of a KvN wave satisfy two equations which 
are completely decoupled and equal to each other. As a consequence, the KvN approach 
to classical statistical mechanics is somehow redundant, since there is not a one-to-one 
correspondence between physical states and waves in the KvN space. As we shall explicitly 
show in the following sections, this redundancy will allow us to calculate classical thermal 
averages in many different ways. 

In order to better understand this point, let us consider the prescription to calculate 
statistical averages within the KvN framework. The mean value of a quantity 0(ip) with 
respect to the probability distribution p(ip, t) can be easily expressed in terms of a generic 
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state \ip a ,t) associated with such a distribution in the following way: 

{0(<p)) p =ty a ,t\0(<p)\il> a ,t) = J dipO(ip)\(ip\4> a ,t)\ 2 = J dipO(ip)p(ip,t), (24) 

where the completeness relation of the eigenstates of (p a has been used. More generally, 
one could also choose to consider a statistical mixture of physically equivalent states \ip a ) 
described by the "statistical operator" 

p(t) = J2Wa,t)(lp a ,t\, (25) 
a 

where the real coefficients f a can be chosen arbitrarily. The mean value of 0((p) with respect 
to p(ip, t) can then be calculated in the following way: 

Thus, every particular choice of the coefficients f a in equation (25) results in a different 
representation of the same statistical distribution p(ip,t). 

Up to this point we have only considered some general aspects of the KvN formulation 
of classical mechanics. Let us now turn to the dynamical aspects and discuss the evolution 
equation in the KvN space, which is the classical analog of the Schrodinger equation. It 
is well known that every probability distribution p((p, t) in phase space evolves in time 
according to the Liouville equation 

id t p{(p, t) = (id q Hdp - id p Hd q ) p(ip, t) = Lp((p, t), (27) 

where H is the Hamiltonian of the physical system under consideration and L is the generator 
of time evolution which is called Liouvillian. The evolution equation for the states ip(ip,t) 
in the KvN space should be consistent with equations (21) and (27), and this requirement 
can be satisfied in a simple way by postulating that ip(ip,t) evolves according to the same 
equation, i.e. 

id t iP(tp,t) = U(tp,t). (28) 

Notice that the KvN waves i^(ip,t) and their square modulus evolve according to the same 
equation because the Liouvillian is a first-order differential operator. This is not the case 
in quantum mechanics, where the wave function and its square modulus satisfy different 
equations because the Hamiltonian is a second-order differential operator. 
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B. Classical Path Integral (CPI) approach 



After this brief introduction to the Hilbert space approach to classical mechanics, we will 
now discuss its path integral counterpart and, in particular, we will derive a path integral 
representation for the kernel 

K fl ^( Vf \e-^-^\ Vl ). (29) 

Instead of following the standard slicing method usually adopted in quantum mechanics, 
we will here proceed in a slightly different and more instructive way. In fact, an implicit 
expression for the kernel can be derived by resorting to the fact that, in classical me- 
chanics, if we follow the evolution of an eigenstate of (p a we get another eigenstate of the 
same operator. 

In order to prove this statement let us first define the derivative operators A = (X q , X p ) 
in such a way that 

(ip\\a\ll>,t) = -id a 4>(<p,t). (30) 

These operators satisfy the algebra 

[<p a ,\ b ]=i6b> [AaA]=0, (31) 

and can be used in order to rewrite the Liouvillian defined in (27) as L = \ a u ab dbH(ip), 
where u ab are the elements of the symplectic matrix [36]. Let us now evolve an eigenstate 
for an infinitesimal amount of time e and then consider the action of (p a on the resulting 
state. By using the commutation rules (22) and (31) we get, up to first-order terms in e, 



ip e \ipi) = e ip \ifi) + 



\<Pi) 



w (p a e~ l Le \ipi)-ie (p a ,L e~ lL£ \ipi) 

= (^ + eu ab d b H( l p i ))e- iLe \ l p i ) «^( £ ;^,0)e-^|^), 

where <f%(t; (p^ 0) is the particular solution of the classical equations of motion (p a = 
uj ab dbH((p) with initial conditions ^ imposed at t — 0. Therefore, e~ tLe \(pi) is still an 
eigenstate of (p a and the associated eigenvalue is <p^i{e; tpi, 0). This result is also valid for a 
finite amount of time and it can be summarized as follows: 

e-^tf-u)^ = \ Vd {t f -^ tl )). (32) 
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As anticipated before, this result can be used to find an implicit expression for the kernel Kfi 
in terms of the solution to the equations of motion. In fact, by exploiting the orthonormality 
of the eigenstates of (p a we get: 

K fi = ((ffle-^f-^ltpi) = ((Pf\(Pd(tf;(pi,ti)) =6((p f -(p d (tf,(p i ,t i )). (33) 

This expression admits a simple physical interpretation, since it basically states that the 
classical probability amplitude of going from ifii to iff in a time interval tf — ti is different 
than zero only if the classical trajectory coming out from ipi at time ti passes through iff at 
time tf. Of course, equation (33) should be regarded as an implicit result, since in most cases 
we are not able to explicitly solve the equations of motion. It is however a good starting 
point for deriving a path integral representation of Kfi. 

To this end, let us first rewrite the kernel as a functional integral of a functional Dirac 
delta involving the whole classical trajectory. This can be done by slicing the time interval 
tf — U in N + 1 subintervals of lenght e and taking the limit N — > oo: 

n r n 

6((pf-(pd(tf,(pi,ti)) = Yl d(p n Y[S((p n+ i- (p cl (t n+1 ;(p n ,t n )) 

n=l J n=0 

^ jWipSbW-MtwuU)). (34) 

In deriving this result we have defined t n = ti + ne and ip n = tp{t n ). Furthermore, V"ip 
indicates that the functional integration is over all the paths ip(t) with fixed end points ipi 
and iff. The functional Dirac delta in equation (34) can be expressed in terms of another 
delta function involving the Hamilton equations of motion: 

<% a - = 5{^ a - uj ab d b H) det{5*d t - u ac d c d b H). (35) 

Both elements on the right-hand side of the above equation can be expressed in terms of 
functional integrals. For the first one we can use the functional Fourier representation 



Sty" - oo ab d b H) = JvXexp^i J f dt X a (v? a - u ab d h H) | , 



(36) 



while the second one, using four Grassmann variables c a = (c q ,c p ) and c a = (c q ,c p ), can be 
expressed as follows: 

det{5^d t - u ac d c d h H) = [ VcVc exp ( - [ 1 dt c a {5^d t - uj ac d c d b H)c b \ . (37) 
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By combining equations (33) through (37) we obtain the following path integral representa- 
tion of Kff. 

K fi = J VipVXVcVc exp ji J ' dt (\ a ip a + ic a c a - ft) j , (38) 
where we have defined 

U = X a co ab d b H + ic a u ab d b d d Hc d . (39) 

The first term in H is just the Liouvillian L, while the second term is responsible for the 
time evolution of the Grassmann variables c a and c a [15]. 

It was shown in [15] that the variables c a turn out to have an interesting physical inter- 
pretation in terms of Jacobi fields, which are the small variations among nearby classical 
trajectories. As a consequence, their correlation functions are related to the Lyaupunov 
exponents of the system described by the Hamiltonian H [37]. Moreover, these variables 
can also be interpreted as a basis of differential 1-forms d(p a [15, 38]. Indeed, notice that, 
due to the anticommuting nature of these variables, the antisymmetry property of the wedge 
product between 1-forms is automatically recovered. Similarly, % itself corresponds to the 
Lie derivative of the Hamiltonian flow [36] associated with the Hamiltonian H. There- 
fore, since the Grassmann variables c a carry some interesting mathematical and physical 
information, it is important to consider the extended Hilbert space [24] associated with the 
path-integral (38). 

From the kinetic term in (38), it is easy to see that the operators c a and c a acting on this 
extended space must obey the following algebra [15]: 

{c a ,h}=5 a b , {e a ,c fe } = {g a ,S b }=0, (40) 

where the brackets above denote anti-commutators. Furthermore, c a and c a commute with 
the operators (p a and A a . In this enlarged Hilbert space, time evolution is generated by the 
operator T-L defined in equation (39), and every state obeys the equation of motion [15] 

idt\^t) = n\^t), (4i) 

which is a generalization of equation (28). The evolution generated by this equation may 
or may not be unitary, depending on how the scalar product is defined. In what follows, we 
will adopt the symplectic scalar product [24] such that (c a Y = ioj ab c b . With this definition, it 
is easy to see that the operator % is Hermitian, and therefore that time evolution is unitary. 
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A brief review of the symplectic scalar product can be found in the Appendix A, while we 
refer the reader to [24] for more details concerning alternative implementations of the scalar 
product. 

C. Universal supersymmetry and superfields 

One striking feature of the approach to classical mechanics we outlined above is the large 
number of auxiliary variables. In fact, while the usual Hamiltonian formulation requires 
2n variables ip a to describe a physical system with n degrees of freedom, the path integral 
approach (or, equivalently, the formulation based on the extended Hilbert space) makes 
use of the 8n variables y? a , A a , c a and c a . However, the redundancy of this description is 
accompanied by a high degree of internal symmetry, as already noticed in [15, 39]. In what 
follows, we will mostly be interested in the symmetry transformations associated with the 
following charges [39]: 

Q H = c a (i\ a + f5d a H{(p)/2) (42) 
On = t a u ab (i\ b - (3d b H(0)/2) , (43) 

where (3 is an arbitrary parameter. By using the algebra given by the equations (22), (31) 
and (40), it is easy to check that these charges commute with the generator of time evolution 
T-L and therefore that they are conserved. Furthermore, by using the Hermiticity conditions 
(c a )t — ico ab c b one can show that Q^ H = iQn, and thus that the operator 

U(e,e) = exp [~eQ H - Q H t) (44) 

generates a unitary transformation (provided e* = it). In order to better understand the 
nature of this symmetry, we can consider the algebra obeyed by the charges Qh and Qh, 
which reads 

{QH,Q H } = -i/3H. (45) 

From this equation, we can see that Qh and Qh are the generators of a non-relativistic 
N=2 supersymmetry (SUSY). Notice that this symmetry has a universal character, in the 
sense that it does not depend on the nature of the physical system, i.e. on the form of the 
Hamiltonian H . 
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This universal supersymmetry allows a simplification of the formalism by combining the 
variables (f a , A a , c a and c a into a single object which bears a close resemblance to the 
superfields one can introduce in ordinary supersymmetric field theories. As it is well known, 
the introduction of a superfield requires us to enlarge the base space. In our case, this means 
that, besides the ordinary time t, two "Grassmannian" partners of time 9 and are needed. 
The triplet r = (t, 9, 9) is called supertime and plays the same role as the superspace in 
supersymmetric field theories. A superfield can then be introduced as follows [23] : 

$ a (t, 9, 9) = y? a (t) + 6c a (t) + 9u ab c b (t) + i99u ab X b (t). (46) 

With our convention for the scalar product 1 , it immediately follows that $ a is Hermitian 
provided 9* = %9. 

Many of the results presented above can be expressed much more elegantly in terms of 
superfields. For example, the algebra given in equations (22), (31) and (40) is fully encoded 
in the much more compact equation [23] 

<l a (t, 9, 9), <S fe (t, 9', 9')] = icu ab 5(9 - 9')5{9 - 9'), (47) 

where we remind the reader that, for a Grassmann variable 9, the Dirac delta is simply 
given by 5(9) = 9. A great simplification occurs also at the level of the path integral. In 
fact, let us consider the Hamiltonian H(ip) and replace the phase space variables ip with the 
superfields $. By expanding H{§) in powers of 9 and 9, we obtain 

H(<S>) = H(ip) + 9N + N9-i99U, (48) 

where the two components iV and N of the supermultiplet above are two functions of ip a , 
c a , c a and A a whose explicit form is not relevant to our analysis. On the other hand, the 
last component of the multiplet is the generator of time evolution defined in equation (39), 
which can be extracted from H(Q) by an integration over the Grassmann partners of time 
9 and 9: 

U = i J d9d9H(<&). (49) 

Similarly, the kinetic term in equation (38) can be rewritten by introducing an integral over 
9 and 9 as follows: 

d9d9 $ p $ 9 + — {X p p + ic p c p ). (50) 



1 We are also defining the complex conjugate of a product of Grassmann numbers as (8162)* = 6\Q\. 
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Because of the surface terms appearing on the LHS of this equation, the path integral (38) 
cannot be expressed solely in terms of superfields. However, these surface terms disappear 
if, instead of considering (which is essentially the transition amplitude between two 
states q) and \ipf, Cf) integrated over the initial and final values of c a ), we calculate the 
transition amplitude between two eigenstates of the superfield $ 9 = q + 8c q + 8c p + i66\ p . 
Notice in fact that all the components of $ 9 commute among each other and therefore they 
can be diagonalized simultaneously (see [23] or Appendix A for more details). In conclusion, 
we have [23] 

= j V"^ q V^ p expji J f dr &>&-H($) j, (51) 

where we have denoted with |<3> 9 ) the eigenstates of & q , and we have defined the differential 
of the supertime r = (t,9,9) as dr = idtdOdO. Notice that, according to our conventions, 
dr is real. Furthermore, we used the notation V"$ q to indicate that we are integrating over 
all the paths $ 9 (t) with fixed end points and 

The formal similarity between the result (51) and the standard path integral expression 
for a transition amplitude in quantum mechanics is striking: 

(qf,tf\qi,U) = J VqVpexp^i jT dt [p? -#(</?)] j . (52) 

It is easy to realize [23] that the "dequantization rules" which let the quantum transition 
amplitude (52) turn into the classical one (51) are the following: 

1) Replace the phase space variables (q,p) with the associated superfields (<3>' 3 , <3> p ); 

2) Extend the time integration to a supertime integration 

dt — > / dr. 



As we will see, in order to establish a connection between classical and quantum field theory 
in thermal equilibrium, we will need to supplement these two rules with a third one, which 
will be derived in section VI B. 



IV. CLASSICAL CTP FORMALISM 



We will now start applying the classical formalism reviewed in the previous section to the 
study of classical thermal field theory (CTFT). In particular, we will first review the path 
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integral formulation of CTFT developed in [14]. As we will see, this approach turns out to 
be related to the high temperature limit of the CTP formulation of quantum thermal field 
theory. 

In this section as well as in the following ones, we will continue to illustrate the general 
formalism by considering the concrete example of a scalar field described by the Hamilto- 
nian 

Of course, our results can be extended to more general systems in a straightforward way. In 
analogy with the previous section, we will indicate the collection of the phase space variables 
as ip — ((p, 7r). 

The starting point of the approach developed in [14] is the following expression for the 
n-point thermal correlation functions [7, 8]: 

... <fr(x n ))p = Z~ x J [difi] U) ... <fr d (x n ; ip h U) e - pH{ipi) . (54) 

In this formula, (f) c i{x; <Pi, U) is the solution of the classical equation of motion associated with 
the Hamiltonian (53), [d(pi] denotes a path integral over all the possible field configurations 
at a fixed time t i; and Zp is the classical canonical partition function defined as 

Zp = J[Me- pHM . (55) 

Equation (54) is just the standard definition of thermal average, and it basically states 
that (4>{xi) ...(f>(x n ))p is given by the product of the classical solutions evaluated at points 
xi, ...,x n with initial conditions ipi averaged over the Boltzmann distribution. 




A. Generating functional 

Equation (54) can be easily re-expressed in terms of an expectation value in the KvN 
space. In fact, the solution of the classical equations of motion can be written as: 



(f) c i(x;ipi,U) = / [d<p f ] (<Pf,tf\<f>(x)\<pi,ti), (56) 
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where x = (i,x) and U <t < tf. This result can be easily proven in a few steps. In fact, 
dropping momentarily the dependence on the space variable x, we have: 

{(Pf,t f \4>(t)\(pi,ti) = J [d(pf][d(p] {tpf,t f \<p,t)<l>{<p,t\<pi,ti) 

[d<pf][dip]6(ip f - <Pd(tf,<p,t)) </>6(<p - <Pci(t;<Pi,ti)) 



[d(p f ] <f) d (t; (p h U) 5 (iff - ipci(tf, ifci(t; tfii, U),t)) 
= <j>d(t;(pi,U), 

where in the first line we have used the completeness of the eigenstates \tp, t) and in passing 
to the second line we have rewritten the transition amplitudes by using equation (33). 
Similarly it is easy to show that the thermal correlation functions introduced in (54) can 
be rewritten as 

... <f>{x n )) p = Z~ l J [d<p f ] [difi] (<p f ,t f \T [0(^0 ... 0(x„)] fa, U) e-WM, (57) 

where ti < < tf for any k — 1, ...,n. Equation (57) establishes a connection between the 
standard approach to the canonical ensemble encoded in equation (54) on one hand, and the 
KvN operatorial approach on the other hand. By resorting to a slicing procedure familiar 
from ordinary quantum field theory and by using the result (38), equation (57) can also be 
expressed as a ratio of path integrals: 

/ Vip VX VcVc 0Qn) ... <j)(x n ) exp {iS - {3H(y (U))} 
(0(x x ) ... <P(x n ))p = JVyVXVcVcexpiiS-PHMU))} ( 5 

where, according to equations (39) and (53), we have that S = Si + S 2 + S 3 with 



/ 



Si = / dt I d 6 x 

■i 

<■*/ 



A^0 + Kit - A^tt + A^ (-V 2 + m 2 ) (59) 
S 2 = I dt I d 3 x [ic<i,& + ic w c n - ic^c* + ic w (-V 2 + m 2 ) c ] (60) 



u 

S 3 j dt / (f'.l 

Jti 



9 \ 1 3 . % 9 - <h 1 2 

31 + — 



(61) 



Notice that we have decomposed S into a part Si, which is quadratic in (ip a , A a ), a part S 2 
which is quadratic in the Grassmann variables (c a ,c a ), and a part S 3 which is quartic in 
the fields and introduces a coupling among the various fields. In particular, the fact that 
is coupled to the auxiliary variables A^, and c n , suggests that a perturbative calculation 
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of thermal correlation functions of the field based on the expression (58) requires the 
introduction of a generating functional Z which depends on four currents, namely 



J VipVXVcVc exp < 


iS-PH{ip{U))+i j 




J 


f VipVXVcVc exp{iS 


-0H(<p(ti))} 



Z = - (62) 



Clearly, the thermal average of a product of n fields can be calculated by taking n deriva- 
tives of Z with respect to and eventually setting the external currents equal to zero: 

l\ n 5 n Z 
i) SJ^Xi) ... 5J^{x n ) 



(0(Xi) ...<p(x n ))p 



(63) 



B. Feynman rules 



We can now use (62) to derive a set of Feynman rules which will allow us to calculate 
thermal correlation functions in a systematic way. The propagators can be easily calculated 
by considering the case in which the coupling g in equation (53) vanishes and, therefore, 
the quartic terms included in £3 are absent. In this case, the Grassmann variables c a , c a 
decouple from <p a and A , and the generating functional (62) reduces to the product of two 
terms, Z = Z 1 [J 4> ,A«] Z 2 [r ]4> ,r ] % with 



/ VipVX exp jiSi - f3H((p(U)) + i j (A"A W + J^) 
Zi[^,A1 = ^ } ~ (64) 



J Vc Vc exp < 


[iS, + J 


f (%c* + c„if)j 


J 


f VcVc exp {1S2} 



Hn^] = - -■ ( 65 ) 



Both these quantities involve only Gaussian integrals and therefore can be calculated exactly. 
In particular, Z\ can be easily calculated by using the result 



(66) 



J ' dtpf J VipVX exp |i«Si + i J (A W A* + J^) 

= JVexp(i I 1 d 4 x J^x^x) + i j S d 4 xd 4 x'J^(x)G R (x - x')A n (x')\ , 
I Ju Jti J 

where <t>o(x) denotes the solution of the free equations of motion with no currents and with 
initial conditions ipi at time tj. More explicitly, we have 

0o(x; <p u U) = j ^ (&(p) cos[E p (t - U)] + ^ ME p (t - e^ x , 
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where again E p = a/p 2 + m 2 , while G R (x) is the retarded Green function, which admits the 
following Fourier representation [28]: 

f d 4 p e ~ip-x 

GR{X) = JWY f-m 2 + l e P ° - (67) 
For a derivation of equation (66), we refer the reader to [40], where an analogous result was 
derived in the simple case of a single harmonic oscillator. At this point, all we need to do 
in order to calculate Z\ is to integrate (66) over the initial conditions ipi weighted by the 
Boltzmann distribution and impose the normalization condition Z\ = 1 when the external 
currents vanish. This is just a Gaussian integral which can be calculated exactly, and the 
final result becomes particularly simple in the limit where U — > — oo and tf — > oo. In this 
case, we simply have 



Zi = exp j-i J d A xd A x' [J^(x)Ap(x - x')J^x') - 2iJ lj> (x)G R {x - x')A n (x')} 



(68) 



where 



d 4 p 2vr 2 2 



A ^ X) = J J2^W\ S(P - m)e ~ iP ' X - (69) 
The calculation of Z 2 is somehow subtle, and the details can be found in the Appendix B, 

but the final result is very simple: 



(70) 



From equations (68) and (70) we can immediately read off the Feynman rules for the prop- 
agators in the momentum space: 

GW = * •=A j8 (p) G^=' • = -iG R (p) 

G^ = - - = -iGU(p) G^=- -=0 (71) 

G c ^ =•■•■•►•■•••= -G R (p) G^ c4 , = ~+ ■■ ■■• = G A (p), 

where we introduced the advanced Green function which, in momentum space, satisfies the 
relation Ga(p) = G R (—p). 

Let us now consider the quartic interactions in (53), which in principle will affect equa- 
tion (62) not only through 1S3 (see equation (61)), but also through the Boltzmann weight 
—/3H(tp(ti)). However, the latter dependence becomes negligible in the limit tj — > —00 
provided [8] we turn off the potential adiabatically by modifying the coupling g as follows: 

g^ge-W, (72) 
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with e a positive infinitesimal parameter. Then, the only contribution will come from iS 3 , 
and the Feynman rules for the vertices in momentum space can be simply deduced from the 
expression (61): 




It is instructive to apply the Feyman rules derived above to calculate the one loop cor- 
rections to the 2-point function (<j)(xi)4>(x2)) p in momentum space, which are represented 
by the following diagrams: 




The contribution of the first two graphs is essentially determined by the following loop: 
V ig f d 4 P A / x ig f d 4 P 2tt 2 2 ig f d 3 p 1 

= W wr Mp) = ^J wrw\ ip ~ m) = ^J Wr& P (75) 

where as usual E 2 = p 2 + m 2 . This contribution diverges linearly, showing that divergent 
corrections appear also in the context of classical (thermal) field theory. Notice however 
that this divergence is milder than the quadratic one encountered in the the case of quantum 
thermal field theory [28]. On the other hand, it is easy to see that the last two diagrams in 
equation (74) cancel each other, since 



where in the last step we have used the fact that a "Grassmannian" loop carries an ex- 
tra minus sign. This cancellation is a natural consequence of the universal supersymmetry 
discussed in section III C. Indeed, similar cancellations appear also at higher orders in 
perturbations theory as shown in [41]. Supersymmetry moreover makes it possible to com- 
pare the classical and quantum perturbation theory via the introduction of super- diagrams. 
Further details will appear in [41]. 

Before concluding this subsection, it is worth commenting on a somewhat technical point. 
The careful reader may think that the two diagrams in equation (76) are both identically 
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zero, since they are proportional to Gr(x = 0). However, as already pointed out by Jeon 
[14], this is not the case unless we choose to regularize the Heaviside step function in such a 
way that 0(t = 0) = 0. The regularization of 0(t) affects in turn the way path integrals are 
discretized. In particular, the choice 0(t = 0) = adopted by Aarts and collaborators [7, 8] 
as well as by other authors [5], corresponds to the so-called prepoint (or Ito) discretization. 

The path integral expression (38) for the classical transition amplitude is based instead 
on the midpoint (or Stratonovich) discretization which implies the regularization #(0) = 
lim e _;.o[#(— e) + 0(e)]/ 2 = 1/2. In fact, a careful calculation of the functional determinant 
(35) in the discretized form shows [42] that the result (37) is correct provided the midpoint 
discretization is adopted. If instead one chooses the prepoint discretization, the functional 
determinant is just a constant and there is no need to introduce the auxiliary Grassmann 
variables [43]. However, in the latter case extra terms depending on ip a and A a may appear in 
the continuum limit 2 and perturbation theory may require different diagrams in the discrete 
and in the continuum limit [44]. Of course this is not always the case and it depends on the 
model under analysis. 

It is important to stress that the sum of all one loop corrections (74) to the two point 
function (0(x 1 )0(a; 2 )) / 3 does not depend on the way 0(t) is regularized, by virtue of equation 
(76). Nevertheless, since the Grassmann variables have a physical interpretation as related to 
the Lyaupunov exponents and a geometrical one (see Section III B) , plus they play a crucial 
role in making the universal supersymmetry manifest (see Section III C), we prefer to adopt 
the midpoint discretization corresponding to the regularization condition 0(0) = 1/2. 

C. Comparison with the CTP formalism at high temperatures 

We shall argue now that the approach to classical thermal field theory we just presented 
[14] can be regarded as the classical counterpart of the quantum CTP formalism. In fact, 
we will show that, after a suitable change of variables, the propagators in the quantum 
CTP formalism are approximated by the propagators of classical fields and A,,- defined in 
equation (71) in the high temperature limit. To this end, let us recall the propagators for 
a quantum scalar field with a quartic self- interaction [28], which are characterized by four 

2 A quick calculation of the functional determinant was also carried out in [15], but that result was not 
accurate enough to keep the discretization into account. 
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propagators, 



Gn(p) = p 2_ m 2 + ie + ^M\pW - m 2 ) (77a) 

G u (p) = 2tt [6»(-p°) + n*(|p°l)] S(p 2 - m 2 ) (77b) 

G 21 (p) = 2tt [0(p°) + n B (|p°|)] 5{p 2 - m 2 ) (77c) 

G 22 (p) = — V— + 2™ B (\p°\)5(p 2 - m 2 ), (77d) 



where % denotes the Bose distribution 

.(A - . 

- 1 ' 



n B (e) = -J—, (78) 



If we now consider the field redefinition 3 



0! \ / 1/2 1/2 \ / 0! 

02 J -1 1 U 



= C/0, (79) 



the propagators transform as G^- = U™ Uj n G mn , and they read: 



Gn(p) = 2tt [1/2 + n B (b°|)] 5(p 2 - m 2 ), (80) 
G 12 (p) = -iG R {p), G 21 {p) = -iG A (p), G 22 (p) = 0. 

By comparing these propagators with the classical ones given in equation (71) we see that 

G 12 ( P ) = G^(p), G 21 (p) = G x ^(p), G 22 (p) = G\ n \ n (p). (81) 

Furthermore, the two sets of propagators become identical in the high temperature limit, 
since Gn(p) reduces to G^p) when j3 — > 0. In fact, 

Gn(p) « 2* Q + 5(p 2 -m 2 )^^ 5(p 2 - m 2 ) = G H {p). (82) 

Thus, equations (81) and (82) show that, as we anticipated, the free propagators for the 
quantum CTP formalism are well approximated in the high temperature limit by the classical 
propagators derived in this section. 



3 Notice that the field redefinition (79) differs from the one commonly adopted in the literature [5, 8, 14]. 
This discrepacy can be traced back to the fact that our variables A a differ by a minus sign from the way 
analogous variables are usually defined in the literature. 
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V. CLASSICAL THERMOFIELD DYNAMICS 



In the previous section we reviewed the approach developed in [14], in which thermal 
averages are calculated starting from the following expression: 

... <t>{x n ))p = Zp 1 J [chpi] (j>d(xi; U) ... 4> c i{x n ] (p h U) e ' mLpl) . (83) 

This way of calculating thermal averages is of course very intuitive, since it is based on 
the idea that the initial conditions of a system in thermal equilibrium must belong to a 
canonical ensemble. However, even though we were able to use the Hilbert space formulation 
of classical mechanics to express equation (83) in terms of the expectation value of a time- 
ordered product of fields, equation (57) is not the natural way to calculate statistical averages 
in the KvN formalism. In fact, the average value of an observable A(tp) with respect to a 
statistical distribution p(ip) is more naturally given by 

{A((p)) = (^\A^m, (84) 

where all the information about the statistical ensemble is contained in the state \ip), which 
must satisfy the relation |(0|</?)| 2 = p((p). 

In this section we will introduce an alternative approach to classical thermal field theory 
in which thermal averages are calculated like in equation (84). The implementation of such 
an approach necessarily involves some degree of arbitrariness since, as we pointed out in 
section III A, a probability distribution p(<p) can be described by an infinite set of states 
in the KvN space. There is however a very powerful criterium which will allow us to 
single out a particular state to describe a canonical ensemble. In fact, we will explicitly 
show that there is only one state in the (extended) KvN space which is invariant under the 
supersymmetry transformations introduced in section III C, and that such a state precisely 
describes a canonical ensemble. 

The idea of expressing thermal averages of observable quantities as expectation values 
of operators acting on a physically redundant Hilbert space is also the starting point of 
an alternative approach to quantum thermal field theory known as Thermofield Dynamics 
[19]. As we will see, the correspondence between Thermofield Dynamics and the classical 
approach we will develop in this section is not only qualitative, but also quantitative. In fact, 
we will derive the classical Feynman rules for a scalar field within this alternative classical 
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approach, and we will show that the propagators for and A^ perfectly match the quantum 
propagators of Thermofield Dynamics in the high temperature limit. 

A. Canonical ensemble and invariance under supersymmetry 

We will now show that the (extended) KvN Hilbert space admits only one supersymmetry- 
invariant state, and that this state is associated with the canonical ensemble. For simplicity, 
we will carry out our analysis for a system with a finite number of degrees of freedom, but 
the result can be extended to field theories in a very straightforward way. Let us start by 
rewriting the supersymmetry charges introduced in section III C: 



where (3 is a real parameter. The request that a state \ip) associated with a probability 
distribution p(<p) be invariant under supersymmetry can be expressed in the following way 



In order to explore the restrictions that these conditions impose on the state let us 
consider the associated wave function in the (ip, c) space [24]: 



c ) = ((-«?)•+, (ic«y +, y#> = m<p) + Mv>) c " + M<p) c p + m<p) (87) 



Notice that, according to the notation introduced in appendix A and due to the Hermiticity 
property (c a )t = iu ab Cb, the bra ((— ic p )*+, (ic q )* +, ip\ is a simultaneous eigenstate of (p a 




(85) 





(86) 



and c a with eigenvalues ip a and c a . In this representation the operators A a and c a can be 
implemented as derivative operators, such that 




d 



(88) 



(89) 



As a consequence the request Qh\iP) = can be rewritten as follows: 



= <(-«?)•+, (id 9 )* 




(90) 



c a (d a + f3d a H(v)/2)4>(^c), 
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If we now use the explicit form (87) for ip(ip,c), we see that equation (90) is equivalent to 
the following two relations: 

(d a + pd a H(<p)/2)M<P) = (91) 
(d p + pd p H(<p)/2)rl> q (<p) = {d q + (5d q H^)/2)^). (92) 

The first equation can be easily solved for any form of the Hamiltonian H, and it has the 
following solution: 

Mv) = K Q e-^)l\ (93) 

where K is for now an arbitrary integration constant. If we now assume that in the Hamil- 
tonian of the system p and q do not couple to each other, it is then reasonable to make the as- 
sumption that ip q and i/j p have a factorized form ip q (<p) = f q (q) 9q{p) and ip p (tp) = f P (q) g P (p)- 
Under these assumption, equation (92) is then equivalent to the following set of equations: 

(d p + pd p H(ip)/2)g q (p) = g p (p) (94) 
(d q + /3d q H(<p)/2)f p (q) = f q (q). (95) 

Let us now turn to the second condition that \ip) must satisfy in order to be invariant 
under supersymmetry, namely QhW) = 0. This condition can be expressed in terms of the 
wave function ip(<f,c) as 

o = ((-u?y+,(u«) m +MQH\ii>) 

= <(-«*)*+, (ic*)* +, V \t a u ab (i\ b - Pd b H^)/2) |V> (96) 
= -^LU ab (d b -pd b H( V )/2)4>( V ,c), 
or, by using equation (87), as 

(d a -/3d a H( V )/2)M^) = (97) 
(dp - (5d p H^)/2)^) = (d q - (5d q H^)/2)^). (98) 

The first equation can be solved analogously to (91) and yields 

Mf)=K 2 e^/ 2 , (99) 

where K 2 is for now another arbitrary constant. Next, by using the factorized form of ip q 
and ip p , we can trade equation (98) for the following two equations: 

(d p -/3d p H(<p)/2)g q (p) = g p (p) (100) 
(d q -/3d q H(<p)/2)f p {q) = f q (q). (101) 
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These two equations, together with equations (94) and (95), indicate that both ip q and tp p 
must vanish. In fact, by subtracting equation (100) from equation (94) we get 

f3d p H( V )g q (p) = =*► g q (p)=0 =*► i/> q ((p) = 0, (102) 

and similarly, from the difference between equations (101) and (95) we derive that 

Pd q H(tp)f p (q) = f p (q) = ^) = 0. (103) 

To summarize, what we have shown so far is that invariance under supersymmetry re- 
quires the state ip(<p, c) to have the following form 

xjj(<p, c) = K Q e-WW 2 + K 2 e m * )/2 c q c p . (104) 

We will now see that the integration constants K e K 2 can be determined by considering 
the probability distribution associated with a state of the form (104). An instructive way to 
do that is to rewrite the expectation value of an observable A(ip) on the state by using 
the completeness relations discussed in the appendix A: 

(A(<p)) = &\A(<p)\iP) = J d<pdc q dc p A(<p)(^\<p,c q -,c p -)((-i( J, y+,(ic q )*+,<p\iP). (105) 

Clearly, in order for equation (105) to agree with the usual definition of average value with 
respect to the probability density p((p), the following relation must hold true: 

J dc q dc p {tlj\<p,c q -,c p -)((-ic p y+,(tc q )*+,<p\tP) = p(cp). (106) 

The integrals on the LHS of this equation can be easily calculated by using the fact that 
(ic p *+, (-ic q *) +, (p\ip) = ip((p, c) is given by (104) and that 

= -K 2 e? HM/2 + Koe-W^ 2 (Tc q *. (107) 

Equation (107) can be easily proved by using the results we derived in appendix A, and can 
be combined with equation (106) to get: 

\ Ko f e -m<P) _ \K 2 \ 2 e^ = p(tp), (108) 
31 



This result shows that, unless the Hamiltonian is bounded from above or /3 = 0, the only way 
to have a probability distribution p(ip) which is positive definite and can be appropriately 
normalized is to have K 2 = and /3 > 0. Then, the value of K is fixed by the request 
p(ip) is normalized to one, and the final result is that the only supersymmetric invariant 
state in the extended Hilbert space which corresponds to a physical probability distribution 
is associated with the canonical ensemble: 

M^c) = Z^ 2 e-^/\ (109) 

where the parameter /3 can now be interpreted as the inverse temperature. Notice that our 
analysis does not depend on the number of degrees of freedom described by the Hamiltonian 
H(ip), and therefore this result is also valid in the case of classical field theories as long as 
there is no coupling between configuration and momentum variables. 

Thus, the fact that essentially the only state which is invariant under supersymmetry 
is associated with the canonical ensemble allows us to interpret the requirement of super- 
symmetry invariance as an algebraic characterization of a system in thermal equilibrium 
[39]. 



B. Generating functional and Feynman rules 

We are now going to use the state (109) to develop an alternative approach to classical 
thermal field theory, in which thermal correlation functions of a scalar field <f> are calculated 
according to the formula 

(0(^0 ... 4>{x n ))p = (4>p\T[4>( Xl ) ... 4>{x n )\ |^>, (no) 

where is the supersymmetric invariant state derived in the previous section. Clearly, 
the correlation functions above can be obtained by deriving the generating functional 

Z= (^Texpji J f dtd 3 x^$ + A n \ n -ifj 4) c' t, -iB n r] 7r ^ |^> (111) 

with respect to the current as in equation (63). As in the case of the classical CTP 
approach, the currents coupled to A*-, and c n are necessary to set up perturbative cal- 
culations. Notice also that, for now, we are considering some finite initial and final times 
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U e tf, but eventually we will take the limits U — > — oo and tf — > +00. By inserting two 
completeness relations at the instants U and tf, we can rewrite Z the following way: 

Z = J[d Vf }[dc f }[d Vl }[dc^^p\ Vf ,c^-,c)-,tf) x (112) 
x (icj*+, (-icj* )+, </?/, t/|T[...] cf-, *i) (-icf)+, ^, t^). 

Each of the three elements in the integrand can be calculated in a straightforward way Let 
us start by noticing that the state is stationary because it is an eigenstate of % with 
eigenvalue 0. In fact, this can be easily shown by using equation (45) to express % in terms 
of the supersymmetry charges: 

%7J> = A {Qh,Qh}\^p) = 0. (H3) 
Then, the first term in the integrand in equation (112) can be easily obtained by taking the 

— 1/2 

complex conjugate of equation (107) with K 2 = and K — Zp , and using the fact that 
does not evolve in time: 

(^\cp f: cj-c}-t f ) = foleMf^c*-,*}-) (114) 

= (^|^,4-,^-) = V 2 e-^)/ 2 c^. 

The last term in the integrand is by definition ij)p((pi, q) and therefore, according to equation 
(109), we have 

«*+, (-i<?)+,<PiMl>e) = Z- Xt2 e~^l\ (115) 

Finally, the second term in the integrand in equation (112) can be calculated using the 
standard slicing procedure together with the path integral expression for the transition 
amplitude (38). Therefore, an opportunely normalized path integral expression for Z reads: 



J Vji exp | 




■[H(<p(U)) + H(<p(t f ))]+i j 




> - ir}^ - ic n rf) j 


J 


f Vfj, exp | 




\ 



(116) 

where V>[i = VfVWcVc, and S is given by the sum of the three terms in equations (59), 
(60) and (61). Notice that the factor c^cj was absorbed in the measure of integration, which 
is now equal to the one used in the previous section. However, the generating functional 
(116) is clearly different from the one of the classical CTP formalism (see equation (62)) 
and as such it will give rise to a different set of Feynman rules. 
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Notice that the only difference between equations (62) and (If 6) consists in the weights 
associated with the initial conditions ipi to the final ones </?/. For this reason, the Feynman 
rules for the propagators of the Grassmann variables will be equal to the ones we derived for 
the classical CTP formalism, see equation (71) and Appendix B. Furthermore, in the limit 
where tj — > — oo and tf — > oo, also the Feynman rules for the vertices remains the same as 
in equation (73), since we are assuming that the interaction is turned off adiabatically in 
the far past and future. Thus, the only Feynman rules that are different with respect to the 
ones derived in the previous section are those for the propagators of and A,,-, which follow 
from the "reduced" generating functional 

! V V VX exp li Si - ^[H(ip(ti)) + H(ip(t f ))] + i [ (A^ + J^)j 

Zi[J*A*M 7-^ f B 1 ■ (U7) 

J V^VX exp - ^[H{ip{U)) + H(<p(t f ))]j 

This generating functional can be calculated by using the result 4 

j VipVX exp + i J (A"A U + J^>) j = (118) 

= J\f6((pf - <p d (tf,(pi,ti))exp\i j d 4 x J<f,(f> c i 



where <p c i are the classical solutions to the equations of motion associated with the modified 
Hamiltonian H = H + j A u [40]. Then, the integration over tpf in (117) becomes trivial, 
while the one over the initial conditions ifii is just a Gaussian integral. The final result is 
the following: 

ZifJ^A 71 "] = exp j-i J d 4 xd 4 x' [J ( j > (x)G (j>( j > (x - x')J^{x') + 

+ 2J 4> (x)G (t>Xn (x - x')A"{x') + A*(x)G^(x - x')A*(x' 
with 

G u = . •=A /3 (p), G^=. . = -l[G R (p) + G A (p)] (119) 

G K4> = . = - l - [G R (p) + G A (p)\ , G x ^ = = A p (p) 



4 Notice that equation (66) can also be derived from (118) by simply integrating over iff. 
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and A/3 denned as in equation (69). These Feynman rules are different from the ones derived 
in the previous section and displayed in equation (71). Notice however that the propagator 
G^ is still the same. Clearly this had to be the case, since G^[x — x') is just the average 
value {4>(x)<j)(x')) p calculated in the free case and, as such, it is a physical quantity that 
cannot depend on the way thermal averages are implemented. On the other hand, there is 
no reason why the Feynman rules for the other propagators should also remain the same, 
since is not an observable quantity. In particular, the most striking difference between 
the Feynman rules in (71) and the ones we just derived is that the propagator Gx^x^ is no 
longer zero. 

C. Equivalence to the classical CTP formalism 

In order to convince ourselves that the approach to classical thermal field theory we just 
developed is physically equivalent to the classical CTP approach, we will now reconsider the 
1-loop correction to the 2-point function ((j)(j))p, given by the sum of the following diagrams: 




+ 



(A) (B) (C) (D) 

If we denote by S / g(m 2 ) the loop made out of G^, which we already calculated in equation 
(75), then we can easily calculate the sum of the first two diagrams by using the Feynman 
rules we just derived: 

(A) + (B) = A p {p)^{m 2 ) (- l -[G A (p) + G fl (p)]) - % -[G A (p) + G R (p)]^(m 2 )A p (p) 

= -iAf,(p)[G R (p) + G A {p)]^{m 2 ). (120) 

This result agrees with the one we would obtain by using the Feynman rules (71) and (73) 
for the classical CTP approach: 

(A) + (B) = Ae(p)^(m 2 )[-iG A (p)] + [-iG R (p)]^(m 2 )Ap(p) 

= -iAe(p)[G R (p) + G A (p)}Zp(m 2 ). (121) 
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Moreover, it is easy to see that the remaining diagrams (C) and (D), cancel each other also 
with this new set of Feynman rule, since the relation (76) still holds true: 



1 <l1 " f -^G A (p) + G R (p)]) (122) 




2 J (2vr) 4 V 2 



4 [J (2tt) 4 nx " J (2tt) 4 2 7 (2vr) 4 

This simple calculation illustrates how the formalism developed in this section yields the 
same physical results as the classical CTP formalism [14], even though the Feynman rules 
are slightly different in the two cases. 



D. Comparison with the Thermofield Dynamics at high temperatures 

We have seen in sections IV C that, at high temperatures, there is a very important con- 
nection between classical and quantum thermal field theories. In fact, after a suitable change 
of variables, the free propagators of the quantum CTP formalism can be approximated by 
the classical propagators (71) for and in the limit j3 — > 0. 

However, the CTP formalism is only one of the possible implementations of quantum 
thermal field theory, as reviewed in Section II. For this reason, one might wonder whether 
the alternative formulation of classical thermal field theory we introduced in this section is 
related to any known quantum formalism. We will now show that this is the case, and that 
the propagators (119) are connected to the propagators of quantum Thermofield Dynamics 
[19] in the limit (3 -» 0. 

Let us start then by summarizing the Feynman rules for the TFD propagators, which 
can be easily obtained from equations (19) with a = 1/2: 

Gn(p) = ~ 2 l ^r— + 2nn B (\p°\)5(p 2 - m 2 ) (123a) 

G 12 (p) = 2tt e^ p ° l/2 n B (\p°\)5(p 2 - m 2 ) (123b) 

G 21 (p) = 2ir e^°\ /2 n B (\p°\)5(p 2 - m 2 ) (123c) 

G 2 2(p) = — 2 ~ 2 - + 2i mB (\p°\)5(p 2 -m 2 ). (123d) 
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Once again, ub stands for the Bose distribution given in equation (78). If we now introduce 
the same change of variables we considered in section IV C, namely 



(124) 



the propagators of the TFD approach transform according to = U i m Uj n G mn into: 

G n (p) = |n B (b |)(e^°'/ 2 + l) 2 5(p 2 -m 2 ) 

G 12 ( P ) = - t -[G R (p) + G A (p)} (125) 

G 21 (p) = -l[G R (p) + G A (p)] 

G 22 (p) = 2nn B (\p Q \)(e^-l) 2 5(p 2 -m 2 ). 

Notice that the mixed propagators are independent of the temperature and are equal to the 
corresponding classical ones : 

G 12 ( P ) = G^(p), G 21 (p) = G x ^(p). (126) 

The other two propagators can instead be linked to the classical ones (119) in the limit 
(3 ->■ 0: 

7T 1 9-7T 

0ll{p) w 2 W\ (1 + 1)2 6{p2 ~ m2) = W\ 6{p2 ~ m2) = Gh{p) 

G 22 (P) - 27r^^y*( P 2 -m 2 ) = ^*(p 2 -m 2 ) = G^(p). 

Thus, we have shown that the Thermofield Dynamics formalism is the quantum counterpart 
to the classical formalism developed in this section, which will therefore be called classical 
Thermofield Dynamics. 

In conclusion, it is interesting to point out that the relation between classical and quan- 
tum Thermofield Dynamics is not only limited to the propagators, and it has a far deeper 
conceptual nature. In fact, quantum Thermofield Dynamics was originally conceived as an 
operatorial approach in which thermal averages of quantum fields are calculated as expec- 
tation values over the so called thermal vacuum, a state belonging to an enlarged Hilbert 
space and such that 

(A) p =(Q,fi\A\Q,fi). (127) 
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Clearly, the ideas behind the quantum Thermofield Dynamics approach bear a close similar- 
ity to the KvN approach to classical mechanics we relied upon throughout this section. From 
this point of view, the existence of a connection between classical and quantum propagators 
like the one we outlined above is not very surprising after all. 

VI. CLASSICAL MATSUBARA FORMALISM 

In the previous two sections we introduced two different approaches to classical thermal 
field theory which, at high temperatures, turned out to be related to two well know real 
time formulations of quantum thermal field theory, namely the CTP formalism and the 
Thermofield Dynamics. This classical-quantum connection was possible because, in both 
cases, the number of variables needed to set up a path integral formulation were twice the 
number of physical degrees of freedom (without considering the Grassmann variables c a 
and c a ). 

However, in section II we reviewed a third approach to quantum thermal field theory, 
namely the Matsubara formalism [18]. This approach is less powerful than the real-time 
formalisms, since it can be used almost exclusively 5 to calculate properties of a system in 
thermal equilibrium, but its formulation does not require the introduction of additional 
degrees of freedom. It is then natural to wonder whether it is possible to introduce a third 
formulation of classical thermal field theory which is connected to the quantum Matsubara 
formalism at high temperatures. As we will see in this section, such a formulation exists 
and shares many formal aspects with its quantum counterpart, namely: 

• it is a "trace formalism" ; 

• thermal averages are implemented via the statistical operator p = e~@ H ^; 

• the statistical operator p can be interpreted as a time evolution operator in a suitably 
defined imaginary direction (more on this later). 

A further and very important feature of this new approach is that it relies heavily on the 
superfield formulation introduced in section III C, and as such 

5 Analytical continuations are possible, although they often involve technical subtleties. 
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• it does not require (formally) a doubling of "variables", since it works directly with 
the supermultiplets $ a instead of dealing with the variables (f a and A a separately. 

For all these reasons, we will name this formalism classical Matsubara approach. 
A. Statistical operator and imaginary time evolution 

As we already pointed out in section III A, the KvN formulation of classical mechanics 
is physically redundant, since every statistical distribution in phase space is associated with 
an infinite set of states in the KvN space. For this reason, although expectation values 
were originally defined as {0{ip)) p = (ip\0\ip) with \ip) one of the infinite number of states 
associated with the probability distribution p(y?), we could also calculate them as follows: 

(OM>p=^P, p = Y,fa\1>a)(rl>a\, (128) 

where the states \ip a ) are all associated with the probability distribution p(<p), and the 
coefficients f a are real. This redundancy in the description of a physical system is ultimately 
responsible for the existence of several equivalent approaches to classical thermal field theory. 
In this section, we will explore an approach based on the Gibbs statistical operator: 

P = e -m<p). (129) 

The fact that this is the same statistical operator used in quantum thermal field theory 
is very intriguing. However, it is important to stress that the properties of the statistical 
operator (129) depend crucially on the Hilbert space it acts upon, and since the Hilbert 
spaces of classical and quantum mechanics are different, classical and quantum correlation 
functions will in general be different even though they are calculated by using the same 
statistical operator. A simple illustration of this fact is provided by the trace of the statistical 
operator (129), which is generally finite for a quantum system, but it is divergent 6 in the 



6 Dealing with statistical operators with divergent trace requires particular attention. For our purposes, 
the most relevant difficulty one encounters is that the cyclicity of the trace may fail. This will not be a 
problem in this section since the approach we will describe is aimed at calculating expectation values of 
fields (p a all at the same time, which commute among each other. However, the fact that the cyclicity of 
the trace may fail would be a major obstacle in trying to extend the classical Matsubara formalism to a 
sort of classical real time formalism. 
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case of a classical system: 

Trp = J dipe-P H ^{ip\ip) ~ 5(0) = oo. (130) 

In the quantum Matsubara approach, the trace of the density operator p = e -P H ^) can 
be written as a path integral associated with an evolution along the imaginary axis of time. 
This is just a consequence of the fact that the operator H(<p), which appears in the density 
operator p, is also the generator of time evolution in quantum mechanics. In the KvN 
approach to classical mechanics, however, the situation is different, since the generator of 
time evolution is no longer the Hamiltonian H, but the Lie derivative of the Hamiltonian 
flow T-L defined in equation (39). Nevertheless, we can still consider the action of the Gibbs 
operator p on the various operators defined on the KvN Hilbert space. In fact, p commutes 
with (p a , c a and c a , but it induces a change in the operators A a , according to the following 
formula: 

\' a = p-^\ a p = \ a + ip u J?-. (131) 

From the previous equation, we can easily evaluate the change induced by the operator p 
on the superfield $ defined in equation (46): 

- d(b a 
p- l <$> a {t,6,6)p = ip a {t) + 6c a (t) + 6u a %{t) + iddu ab \ b {t) - 



dt 

= fi a (t - (566) + 6c a (t) + 6u a %(t) + iddco ab \ b (t) 

= &{t- $66,6,6). (132) 

According to our conventions [23] about the Grassmann partners of time, 6 and 6, the 
product 66 is purely imaginary. Therefore, we can think of "complexifying" the time variable 
via the following rule: 

* ->■ t 1 + 66t 2 , (133) 

where, for the purposes of this paper, it will be sufficient to take t\ and t 2 as real variables 7 . 
The action of p on a superfield can then be interpreted, according to equation (132), as a 
"time" evolution along the 66 component in the (ti,t 2 ) plane, as shown in Figure 3. 

It is also interesting to consider the action of p on states in the KvN Hilbert space, and in 
particular on the eigenstates of $^ in the Heisenberg picture. If we multiply the Hermitian 



7 Notice however that there is a fundamental difference between the usual imaginary unit i and the product 
66, namely that i 2 = -1 while {66f = 0. 
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t-2 



FIG. 3: Time evolution along the 99 direction. 

conjugate of the eigenvalue equation 

t\ ¥(t, 6, 6) = t\ $*(0, 6) (134) 

by p and we insert the identity pp^ 1 = 1 after each bra, we easily get 

t\ p &(t - (366, 6, 6) = t\ p $ (0, 6). (135) 

This equation clearly implies that t| p is an eigenstate of the operator $^(i — /36>0, 0, 0), 
which can be indicated as follows: 

($*,t|p = ($*t-£00|. (136) 

The analogy with the quantum Matsubara formalism is striking, and equations (132) and 
(136) are the main results of this subsection which will allow us to set up a path integral 
formulation based on the statistical operator p = e~^ H ^ . 



B. Path integral 

We now have all the necessary ingredients to develop a path integral formulation for the 
classical Matsubara approach. The first step is then to calculate the trace of p on the basis 
of the eigenstates of ^ as shown in equation (A23) of appendix A. By virtue of equation 
(136), we get: 

Trp = J[d&+\ ($ ,O|p|<i^,O) = y*[d$*]($* £00|$* 0), (137) 
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where we have introduced a new superfield which differs from the usual one <3>^ defined 
in (46) only for the sign in front of the Grassmann variables and c w : 



¥(x, 9, 9) = (f)(x) - 9c*(x) - 9c n (x) + i99\ 7T (x) = $*(ar, -9, -9). 



(138) 



Notice also that, for later convenience, we have inverted the order of 9 and 9 in (137) 
compared to equation (136). By following the standard slicing procedure, equation (137) 
can then be rewritten in terms of a path integral. To this end, let us divide the interval (399 
in N + 1 subintervals of length e99 and insert N times the completeness relation (A31) for 
the eigenstates of The transition amplitude in (137) becomes then: 

N N 

($+,poe\$+,o) = n / [d<] Yl($i + i,sn + iemis n ee), (139) 



n=l 



n=0 



where s = 0, sn+i = /3, s n +i — s n = e and we have chosen $q — ^> &n+i = According 
to the definition (46) of superfields, we can rewrite the Hamiltonian H(<p) in terms of the 
superfields as: 



H 



m - j 



id9d9 



i99H($) =H($) 



(140) 



where the notation | indicates the first component of the multiplet H(&), which can be 
also obtained by setting 9,9 = 0. By inserting a completeness relation for the eigenstates of 
in each of the infinitesimal amplitudes appearing in (139) we obtain: 



(<S>i +1 ,s n+1 99\t>is n 99) = (Oe-^W 

« (^+i|exp [-6 j d 3 x 



7T\2 



exp 



-eV(*%) 



[d®l\ ex P 



, (Ml) 



where we have introduced the discretized time derivative $^ = (<&^ +1 — as well as the 

Hamiltonian density Jif for the fields associated with the Hamiltonian H . By combining 
equations (137), (139) and (141) and turning to the continuum limit, we finally obtain the 
following path integral expression for the trace of the statistical operator p: 









*J! ds J 







(142) 
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where the path integral is over trajectories in the space of the superfields which satisfy the 
boundary condition $^(/%>0, x, 0, 9) = $^(0, x, 6>, 9). The result (142) can also be rewritten 
in a simpler way by observing that the integral over the parameter s can be converted into 
an integral over the time variable t(s) = 99s. Using the standard rules of calculus, we get: 

d& 



99 



99ds = dt. 



ds dt ' 

The final result for the trace of the density operator is the following: 

.d^ 



Tr p = / V&V®* exp 



J drj d 3 x U 



dt 



(143) 



where dr = idtd9d9. This is the path integral for the Matsubara approach to classical 
thermal field theory. It is interesting to compare this expression with the path integral for 
the quantum Matsubara approach, i.e. 

-ip 



Trp 



V(pV7i exp 



r-ip r 

7. dt J' 



(144) 



The formal analogy between (143) and (144) is striking, and in fact the quantum expression 
can be turned into the classical one by using the two dequantization rules [23] already 
mentioned at the end of section III C, together with a new rule: 



3. i(3 



99(5. 



C. Generating functional and Feynman rules 



We can now take the path integral expression (143) as our starting point to systematically 
derive the classical perturbative calculations in a way which closely resembles what is usually 
done in the quantum Matsubara formalism. In fact, we can define the generating functional 
as follows: 



Z[J*] 





''J 


-pee / d ^<t> \" 
da \$* d ($) + J^j 




r m ( d& \ 





(145) 



where we have simplified the notation by defining a = (x, 9, 9) and da = i d 4 x d9 d9, and by 
introducing the following supercurrent: 



J (a) = -A n (x) - i9rf{x) + i9fj 4> (x) - i99J^(x). 
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(146) 



As usual, expectation values can be calculated systematically by taking functional derivatives 
of Z[J^] and setting = at the end. However, it is important to realize that, this 
procedure allows us to calculate only time-independent expectation values - another point 
of contact with the quantum Matsubara formalism. In fact, if we parametrize the path in the 
"complexified" time plane shown in Figure 3 as t(s) = —s66, then the functional derivatives 
w.r.t. a supercurrent evaluated on this path can be naturally defined in such a way that 



(147) 



5J^t(s'),^,e',e') 

Then, by using this equation together with the Taylor expansion $^(t(s), x, 9, 9) = 
$^(0,x, 9,9) — sOO <E>*(0, x, 0, 0), we get the following result for the first derivative of the 
generating functional evaluated at = 0: 



1 5Z 



i 530(a) 



5J<f,(a) 



da'J^a'^ia')) = 99 (0(0, x))^ (148) 



This result clearly shows that every functional derivative w.r.t. J </,(<%) will contribute with 
a factor of 0(0, x). Therefore, by using the generating functional (145), we will only be able 
to calculate thermal correlation functions at time t — 0. In fact, by taking n derivatives 
of the generating functional w.r.t. the supercurrents and setting = 0, we can get the 
expectation value of the product of n fields <fi evaluated at different space points, but at the 
same time t — 0: 



iSJ^otx) idJ^atn) 



\\9i9~i 



i=i 



(0(O, Xl )---0(O,x n ))/ 



(149) 



A good way to see that these are indeed the only quantities we can calculate by using 
the generating functional (145) is to realize that the coupling term between the external 
supercurrent and the superfield $^ actually reduces to 



da 3^ 



rf 3 xA^(O,x)0(O,x) 



(150) 



We can now use the generating functional (145) to derive a set of Feynman rules for a scalar 
field with quartic self-interaction described by the Hamiltonian (53). As we will see, these 
rules will involve only one propagator and one vertex, exactly like in the quantum Matsubara 
formalism. 
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As usual, the propagator can be determined by neglecting the quartic interaction. To this 
end, it is more convenient to work in the operatorial formalism where, up to a normalization 
constant, the free generating functional reduces to: 



Z [J«d ~ Tr 


( r" 

pTexp % I 


ee \ ~ 
da 









J [dip] exp | -(3 J d 3 x 



Notice that, in passing from the first to the second line we have used equation (150). The 
integral over n in equation (151) yields an overall multiplicative constant, while integral 
over is a standard Gaussian integral that can be calculated by considering the Fourier 
transform of the field <fi and the current A 71 " w.r.t. the spatial variables only. By imposing 
the normalization condition Z [J</, = 0] = 1 one finally gets: 



ZoW = exp{-^ J ^A*(0,p)^[A*(0,p)]*} 



i r pSe r d 3 v i 

< ^>\-2j dTdT ' j (2^ J ^ r ' p )^ J ^ r/ ' p )^ ( m 



where E p = \/p 2 + m 2 . Equation (152) clearly shows that the classical Matsubara formalism 
admits only one propagator which is given by 

1 



G^t, p ) = — + 0(t). (153) 



Although equation (152) clearly shows that terms of order 0(r) or higher in the propagator 
Gp(r,p) do not contribute to classical thermal correlation functions, we have included such 
terms in (153) because they play a role in the comparison between G c Jt, p) and the propa- 
gator of the quantum Matsubara formalism at high temperatures. In fact, according to the 
general result (12), the propagator for the quantum Matsubara formalism in the (t, p) space 
and in the high temperature limit takes the form 

G,(t, P) = =§gi («*' + e<^--»<) ^ ^ = J_ + 0(f ). (154) 

Hence, we see that the free quantum propagator correctly reduces to the classical propagator 
in equation (153) in the high temperature limit. 

For completeness, let us mention that the rule for the vertex can be easily derived by 
using the fact that the full generating functional can be expressed in terms of the free one 

45 



as 

-pee 

Z[J ] = exp 



-i 



drd 3 xV / idd'dO 1 




x , , \ . Z„[J ]. (155) 

where V is the perturbation, which in general contains terms of order higher than quadratic. 
In particular, in the case of a quartic interaction with coupling g we obtain the following 
rule for the vertex in the (r, p) space: 

/ p-pGO 

~ig / dr. (156) 
Jo 

D. Equivalence to the classical Thermofield Dynamics and CTP formalism 

In conclusion, it is interesting to check that the classical Matsubara formalism we de- 
veloped in this section reproduces the same static results that can be obtained from the 
classical CTP formalism or Thermofield Dynamics. To this end, let us consider once again 
the one loop corrections to the two-point correlation function. Within the classical Matsub- 
ara formalism, these corrections are described by one single diagram, and according to the 
rules (153) and (156) in the (t, p) space they are equal to 



' r ° G}(n-r,k)G}(0,p)G^r-r 2 ,k) 




2(3 2 E* J (2tt) 3 £2- ( 157 ) 

This results should be compared with the one we already calculated in section VC and 
shown in equation (121). By taking the Fourier transform w.r.t. 0-th component of the four 
momentum k and setting the time equal to zero we get 

W + (*) = / i~ iA ^ G ^ + | / (gi ( 158 ) 

and thus we recover the same result we obtained within the classical Matsubara formalism. 
This very simple example is sufficient to illustrate how, in general, static results can be 
obtained much more easily within the classical Matsubara formalism, since fewer diagrams 
are needed. In our example, the single diagram (157) should be contrasted with the four 
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diagrams shown, for example, in equation (74). Nevertheless, the classical Matsubara for- 
malism has its limitations, and when it comes to calculating dynamical quantities, one is 
forced to adopt one of the classical approaches discussed in the previous sections. 

VII. CONCLUSIONS 

In this paper we studied the connection between classical and quantum thermal field 
theory at high temperatures. To this end, we worked, at the classical level, within the 
framework of the KvN and CPI formulations of classical mechanics. 

After reviewing briefly the most common formulations of quantum thermal field theory 
and the KvN and CPI approaches to classical mechanics, we presented the path integral 
approach to classical thermal field theory explored in ref. [14]. We explicitly derived a 
set of Feynman rules to systematically calculate classical thermal correlation functions of a 
classical scalar field with quartic self-interaction, and we showed how the free propagators 
of the quantum CTP formalism are well approximated by the classical propagators at high 
temperatures. We also improved upon the work done in [14] by showing how the cancellations 
occurring among classical diagrams can be seen as a natural consequence of a universal N = 2 
supersymmetry. 

In the second part of the paper, we developed two alternative approaches to classical 
thermal field which turned out to the related to the quantum Thermofield Dynamics and 
Matsubara formalism at high temperatures. In both cases, we derived the Feynman rules 
for a classical scalar field with quartic self-interaction, and we explicitly showed how the free 
quantum propagators reduce to the classical ones in the high temperature limit. 

The equivalence of these three approaches to classical thermal field theory was illustrated 
by considering the one loop corrections to the 2-point thermal correlation function. A more 
in depth study of higher order corrections and the cancellations which occur among classical 
diagrams due to the universal supersymmetry will be the main focus of a follow up paper 
to appear shortly [41]. In that paper, it will be shown how the use of super- diagrams in the 
classical Matsubara formalism greatly simplifies classical calculations. 

The fact that the classical Matsubara formalism is the only one that is entirely formulated 
in terms of superfields is actually puzzling, and the development of a superfield formulation 
for the classical TFD and CTP approaches would further strengthen the formal analogy 
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between classical and quantum thermal field theories. This is however a topic which we 
leave for future study. 
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Appendix A: The Symplectic Scalar Product 

We will now briefly review the definition and some basic results concerning the symplectic 
scalar product. For more details, we refer the reader to [24]. 

Let us start by introducing some basic formulae which do not depend on the way the 
scalar product is implemented. For simplicity, rather than considering the full extended 
Hilbert space of classical mechanics, we will restrict ourselves to the subspace spanned by 
c a and c a . Because of the algebra (40) satisfied by such operators, there must exist a state 
|0— , 0—) such that 

c 9 |0-,0-) = c"|0- 0-> = 0. (Al) 
By letting the operators c a act on this state, we can then define three more states as follows: 

|0+,0-) = g,|0-,0-) (A2) 
|0-,0+) = — cp|0— ,0— > (A3) 
|0+,0+) = g p |0+,0-) = c p c 9 |0-,0-). (A4) 

The four states introduced above form a basis for the Hilbert space on which c a and c a 
are acting. Equation (Al) shows that the state |0— ,0— ) can be regarded as a simultaneous 
eigenstate of c q and c p with zero eigenvalue, but we can also consider more generic eigenstates 
of these two operators with Grassmann eigenvalues a q and a p : 

\a q -,a p -) = e^ 9 ^e" aPSp |0-,0-). (A5) 
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Similarly, we can also consider eigenstates of the operators c a . For example, the state 

\a q -,a p +) =e- aqt «e- a *> dP \0-,0+) (A6) 

is a simultaneous eigenstate of c q and c p . 

Let us now consider the Hermitian conjugate of the expressions (Al) through (A4) and 
denote with (0 + , — | the bra dual to the ket |0 + , — ) (and similarly for the other bras). By 
using the Hermiticity properties of c a and c a with respect to the symplectic scalar product, 
namely (c a )^ = iu ab c b , we easily obtain: 

(0-,0- |c p = (0-0- \c q = (A7) 

(0+,0- I = -i(0-,0- \c p (A8) 

(0-,0+ | = -i(0-,0- \c q (A9) 

(0 +, + | = i (0 +, - \c q = (0 -, - \c p c q . (A10) 

If we now impose the normalization condition 

(0+,0+ |0+,0+) = 1, (All) 

we can calculate all the scalar products among the states which make up the basis of the 
Hilbert space by solely using the algebra (40) of the operators c a and c a . It turns out that 
the only non-vanishing products are 

(0-,0 + |0+,0-) =i (A12) 
(0+,0- |0-,0+) = -i (A13) 
(0-0- jO — , — > = -1. (A14) 

Equation (A14) clearly shows that the symplectic scalar product is not positive definite. 
This is unfortunately an unavoidable feature of the (extended) KvN formalism: in fact, it 
was shown in [24] that the requirement that time evolution is unitary necessarily implies a 
scalar product which is not positive definite, and vice versa. Nevertheless, this result does 
not cause any trouble at the physical level and it has an interesting physical interpretation 
[24]. Because of equations (All) through (A14), the completeness relation reads 

|0+,0+)(0+,0+ I -i|0+,0-)(0-,0+ I + 

+i|0-,0+)(0+,0- | - |0 — , — )<0 — , — I = 1 (A15) 

49 



and, consequently, the trace of a Grassmann-even operator A is 

Tri = (0+,0+|i|0+,0+)-i(0-,0 + |i|0+,0-)+i(0+,0-|i|0-,0+)-(0-,0-|i|0-,0-). 

(A16) 

By using equations (All) through (A14), it is also possible to calculate the scalar product 
between eigenstates of c a and c a . In particular, three useful results that we used in this paper 
are [24] 

((*&)*-, (iP q )* + \a q -, a p +) = 1 5{(5 q - a q ) 5((3 P - a p ) (A17) 
(H/3 P )*+, (-i/3 q )* - \a q -, a p +) = -i exp (f3 p a p + (5 q a q ) (A18) 
((i/3 p )*-, (-i/3 q )* - a p -} = - exp (f3 p a p + I3 q a q ) , (A19) 

where we denoted with (a q —,a p + | the Hermitian conjugate of the ket \a q — ,a p +), and so 
on. Notice that, because of the relation (c a y = iu ab Cb, the bras in equations (A17), (A18) 
and (A19) are eigenstates of (c q , c p ), (c p , c q ) and (c qi c p ) respectively, and the various factors 
of i were introduced, together with the complex conjugate operation, in order to simplify 
the form of the eigenvalues. It is also possible to derive the following completeness relations 
involving the eigenstates of c a and c a introduced above: 

J da q da p \ V , a q -, a p -) ((-ia p )*+, (ia q )* +,<p\ = 1 (A20) 

% J da q da p \a q -, a p +) ((ia p )* -, (ia q )* + \ = 1, (A21) 

and to use these relations to calculate the trace of a Grassmann-even operator A. For 
example, by using equation (A21) we get 

= i J da q da p ((ia p )* -,(ia q Y + \A\(-a q )-,(-a p )+). (A23) 

Notice that, in going from the first to the second line, the sign of the eigenvalues within the 
ket has changed due to the fact that we have interchanged the position of the bra and the 
ket. 

To conclude this appendix, let us apply the general results introduced above to the 
extended Hilbert space of classical mechanics. In particular, it is interesting to see how such 
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results take a much simpler form when they are expressed in terms of superfields [23]. In 
fact, according to our notation and to the definition of the superfields given in equation (46), 
their eigenstates are 

|$*) = 10^,^-,^+} (A24) 
ee |A^,7r,^+,c w ->. (A25) 

Then, the eigen-bras of the superfields are defined as follows 

($*| ee ((f), X n , (icn)*-, {ic' 1 ')* + | (A26) 
ee (A ,tt, (-icT+, (-it-*)* - I, (A27) 

and the scalar products between such eigenstates can be easily deduced from equations 
(A17) and (A18): 

(Sfl**) = (A28) 
= A/"exp ^ <Pxd6d0$*$*^ , (A29) 

where A/" is an irrelevant overall constant and we have introduced the functional delta 

8($t - $t) =15^ - faMKi - X w2 )5(cf - 4)5(c nl - c w2 ). (A30) 

Equations (A21) and (A22) imply that the eigenstates of and must satisfy the fol- 
lowing completeness relations: 

J = J [d$ n ]\$ n )($ n \ = 1, (A31) 

where the integration measures are defined as [d®*] = i[d(f>][dX n ][dc! l '][dc n ] and 
[d^] ee ifrfyrjffiA^jffic^j^c 71 ], and the square brackets indicate that the integrals are car- 
ried out over field configurations at fixed time. Finally, according to equation (A23) the 
trace of a Grassmann-even operator acting on the KvN extended space can be calculated on 
the basis of the eigenstates of ^ as 

Tri = y"[^](^|i|^), (A32) 

where we have introduced a new superfield which differs from the usual one $^ defined 
in (46) only for the sign in front of the Grassmann variables and c n , as shown in equation 
(138). Equations (A28), (A29), (A31) and (A32) are the main results that were used in 
section VI to give a path integral formulation of the classical Matsubara formalism. 
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Appendix B: Free Propagator for the Grassmann variables c a and c a 

In this appendix we will derive the Feynman rules for the propagators of c a and c a by 
explicitly calculating the path-integral: 

z 2 = Mjmvc ex p {« / A ^ - - mv> - -»V - ^ - i^l) , (Bl) 

where Af is determined by the requirement that Z 2 — 1 when r\ a = fj a = 0. As a first step, 
we can rewrite the argument of the exponential in a more compact form as follows: 

% J d 4 xd 4 x'6(x' - x) [c a (x')(i5 a b d t - M\)c\x) + c a (x')P a 6 ^(x) + fj a (x')Q\c b (x)} = 

= i[cDc + cPr] + f)Qc], (B2) 

where we have defined: 

I i\ ( \ 

yi(v 2 -m 2 ) oy \° v 

and 

x) = 5(x' - x) {i5ld t - M\) 

P(x',x) = 5(x'-x)P a b (B3) 

Q(x',x) = S(x'-x)Q\. 

According to the usual Feynman prescription, let us modify the differential operator D by 
adding an infinitesimal parameter e in such a way that D — > D + ie. Equation (B2) can now 
be rewritten as 

i[c + f] Q (D + ie)- 1 ]^ + ie)[c + {D + ie)- 1 P r]] -ifjQ(D + id)' 1 Prj. (B4) 

After introducing new integration variables d = c+fjQ (D + ie)" 1 and d = c+(D+id)~ 1 Prj, 
the calculation of the path integral becomes straightforward, and by requiring that Z 2 be 
normalized to one, we obtain the final result: 

Z 2 = exp(-if ] Q(D + ie)- 1 Pr]) . (B5) 

If we now denote with G the inverse of the differential operator (D + id), the relation 
(D + id) G = 1 can be written more explicitly in Fourier space as follows 

p° + ie -i \ ( 1 , 

2 n G(p) = . (B6) 

%E% p° + tej yo 1 
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This last equation can be put in diagonal form via a change of basis done via the matrix: 

W{p) = I t/Ep " 1/2 I , (detW(p) = 1) . 

In fact, defining Gr>(p) = IU _1 (p)G(p)IU (p), the (B6) becomes: 

p ° +k - e ° W( 10 

p° + ie + E p J yo 1 

From this we can get the solution: 

GmJ^^-^ ° V (B7) 

^ l/|p + it + Bp] / 

The Feynman rules in momentum space are obtained by calculating: 

/ -l/l(p° + ie) 2 -El}\ 

-iQ{p) G{p) P(p) = -iQ(p) W{p) G D (p) W-\p) P(p) = 

V° J 

Notice that the only non-zero compoment is the one that connect the currents rf e fj q . 
Expanding up to the first order in e the term (p° -Me) 2 and re-defining e — > e/2, we get that 
this component is the opposite of the retarted propagator in momentum space: 

G R (p) ' 



p 2 — m 2 + iep" 
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